#include <stdio.h>
#include <stdlib.h>
#include <math.h>
//#include <atlas/clapack.h>
#include <atlas/cblas.h>
#include "test_dgbmv.h"




void testfunc(double matrix[],  double matrix2[], double vector[], double loesung[]){


int m  = 4;
int n  = 4;
int kl = 1;
int ku = 1;

cblas_dgemv(CblasRowMajor,CblasNoTrans,m,n,1.0,matrix,m,vector,1,0.0,loesung,1);
int i;
/*for(i = 0;i<4;i++){
    printf("   %f \n",loesung[i]);
}*/
loesung[0] =0.0;
loesung[1] =0.0;
loesung[2] =0.0;
loesung[3] =0.0;
cblas_dgbmv(CblasRowMajor,CblasNoTrans, m , n,kl,ku,1.0,matrix2,3,vector,1,0.0,loesung,1);
for(i = 0;i<4;i++){
    printf(" hallochen  %f \n",loesung[i]);
}
double s;
s = cblas_ddot(4,vector,1,loesung,1);
printf("loesung %f \n",s);
}
